Gamma Rays from the Tycho Supernova Remnant: Leptonic or Hadronic ? 
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Recent Fermi and VERITAS observations of the prototypical Type la supernova remnant (SNR) 
Tycho have discovered 7 rays with energies E in the range 0.4 GeV E ^IQ TeV. Crucial for the 
theory of Galactic cosmic-ray origin is whether the 7 rays from SNRs are produced by accelerated 
hadrons (protons and ions), or by relativistic electrons. Here we show that the broadband radiation 
spectrum of Tycho can be explained within the framework of a two-zone leptonic model, which 
is likely to apply to every SNR. A model with hadrons can also fit the radiation spectrum. The 
hadronic origin of 7-rays can be confirmed by Fermi spectral measurements of Tycho and other 
SNRs at < 300 MeV. 

PACS numbers: 95.85.Pw, 95.85.Ry, 98.38.Am, 98.38.Mz 



I. INTRODUCTION 

Tycho (G1204-1.4, 3C 10) is the bright spherical super- 
nova remnant (SNR) of SN 1572 [l| with a diameter of 
K, 8'. The distance d to the source is not well-established, 
but is in the range d~ 1.5-4 kpc (e.g., 01 )• The narrow 
^rim ~ 4" spherical rim of hard X-ray emission is coin- 
cident with the outer edge of the radio emission [3] and 
is explained by synchrotron emission of relativistic elec- 
trons accelerated to TeV energies by a forward shock 
moving with speed Vsh ~ 4600 (£i/2.3kpc) kms~^ In- 
frared observations from AKARI Q between 9 and 160 
jjL reveal thermal dust emission in the NE and NW shells. 
The soft X-rays also reveal a rich line structure implying 
significant thermal radiation [Tj. 

Gamma-ray fluxes from Tycho between w 1 and 10 
TeV have been recently reported by the VERITAS col- 
laboration Q . The Fermi L AT Collaboration announced 
a detection in the energy range from « 0.4 to 60 GeV 0, 
with energy spectral index a ~ 1.3. The total Fermi and 
VERITAS spectra are described by a single power law 
with a « 1.1-1.2 from w 500 MeV to 10 TeV. Nonther- 
mal emission at lower frequencies, from radio through X- 
rays, is universally attributed to synchrotron radiation, 
which is a signiture of primary electron acceleration. 

In leptonic models, the 7 rays could result from elec- 
tron bremsstrahlung and/or Compton-scattering pro- 
cesses. In hadronic models the 7 rays arise from the 
decay of, predominantly, 7r° mesons produced in inter- 
actions of cosmic-ray protons with gas and dust in the 
remnant. By modeling its broadband radiation fluxes, 
Morlino & Caprioli |10.] argue that the reported 7-ray 
fluxes provide strong evidence for hadronic acceleration 
in Tycho as the leptonic models seem to fail. A similar 
conclusion was derived by Volk et al. ll| on the basis 
of the then-existing upper limits of the GeV- TeV 7-rays. 
Because of the fundamental importance to the theory of 
SNR origin of Galactic cosmic rays, here we undertake a 
detailed analysis of the robustness of those claims. 

Implicit in the previous modeling is that the detected 



synchrotron and Compton emissions are both produced 
by the very same population of relativistic electrons. 
This is the working assumption of the commonly used 
single-zone model approach. In a SNR environment with 
spatially non-uniform magnetic fields a more realistic 
description of the nonthermal emission should allow at 
least two zones [l^, [3] , where zone 1 represents a strong 
magnetic-field region in the shock vicinity where acceler- 
ation is strong, and zone 2 adjacent (further downstream 
) to zone 1, with larger volume but lower magnetic field, 
into which the zone 1 electrons escape by convection with 
plasma and/or by diffusion. Because the synchrotron 
emissivity in zone 2 is significantly lower than in zone 
1, while the target photon field in both zones is basi- 
cally the same, the Compton 7-ray fluxes in the two-zone 
model can greatly exceed the predictions of the single- 
zone model. 

Within this framework, we show that a leptonic model 
alone can explain well the reported radio through 7-ray 
data of Tycho. The TeV 7 rays are made primarily by 
Compton emission, and bremsstrahlung produces most 
of the radiation at GeV energies. The two-zone model is 
described in Section 2, model parameters and interpre- 
tation of the Tycho's broadband spectrum are given in 
Section 3. As concluded in Section 4, the most impor- 
tant conflrmation of the hadronic cosmic ray hypothesis 
remains the detection of the 7r°-decay feature. 



II. TWO-ZONE MODEL FOR SNRS 

It is widely believed that cosmic rays are accelerated 
through Fermi processes by SNR shocks. Both first-order 
(shock) and second-order (stochastic) processes rely on 
differences of magnetic field in different regions over 
which particle diffuse. Strong magnetic fields favored 
for efficient acceleration can be produced in the vicinity 
of the strong forward shock through non-linear amplifi- 
cation of the magnetic turbulence/ Alfven waves by the 
cosmic-ray ions accelerated at the shock (e.g., [l3. Il5|). 
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The narrow 6*11111 ~ 4" spherical rim of hard X-ray 
emission Q is evidence for such processes operating 
at the forward shock of Tycho. Interpretation of the 
rim width h « 1.93 10~^dkpc pc, at the source distance 
dkpc = d/1 kpc, as due to fast synchrotron coohng of 
^ lOTeV electrons implies a high magnetic field in the 
rim, B > 400 /iG fTTI|. In effect, this explanation leads to 
a single-zone model where all the TeV 7-rays and the X- 
rays are assumed to be produced by the same electrons. 

However, the rim could also be explained as the re- 
gion with strongly enhanced magnetic field behind the 
shock. The rim width then refiects the length scale of 
damping the rtiagnetic turbulence rather thanjjpoling of 
TeV electrons [1^ (see also the discussion in [17]). Rela- 
tivistic electrons accelerated at the rim escape from zone 
1 (the rim) into the zone 2 by diffusion and/or convec- 
tion with the fiuid. The theory of cosmic-ray transport 
is highly developed, and involves at least the spectrum 
and anisotropic character of the turbulence [e.g.,il8., .19.]. 
Nevertheless, reasonably accurate general estimates can 
be made. For Tycho the convective escape time Tc = 
h/v'fi w 40 yr, where ~ Wsh/4 ~ 500 dkpc kms^^ is the 
fluid speed in the forward shock frame. The diffusive es- 
cape time from zone 1 is t^h ~ h'^ /{2k) where k = Ascc/3 
is the diffusion coefficient. The mean scattering length 
Asc of relativistic particles is equal to the gyroradius Vgy^ 
in the Bohm diffusion limit, which is attained if the ratio 
of the turbulent magnetic field fluctuations to the mean 
field is 77 = \SB\/B ~ 1. A value 77 < 1 is probably more 
realistic in general. In this case Age ~ Tgyc/'rf (e.g., fisf). 

The diffusive escape time of electrons with Lorenz fac- 
tor 7 from the rim can be estimated as Tjif ~ 1.1 x 
lO'* 77^ B^G^kpc y^- Comparing this with the syn- 
chrotron cooling time isyn ~ 2.45 x 10^"^ ^^g7^^ one 
finds that accelerated electrons (of any energy) could 
cool in the rim before escaping downstream only if i? > 
1.31 Tj^^/^mG. Even in the limit of 77 = 1 this field is high. 
For any smaller magnetic field, electron escape from the 
acceleration region before cooling cannot be prevented. 
We note that this very effect explains the detection of 
non-thermal X-rays not only from the rim but also from 
a wider region interior to the rim. 

Weaker magnetic fields interior to the thin rim could 
be the result of dissipation of the magnetic turbulence in 
the thermal plasma downstream of the forward shock on 
timescales Tc yr. Therefore the region of the shell interior 
to the rim forms zone 2. This zone extends at least down 
to the contact discontinuity at rco/^'sh ~ 0.927 ^4], in 
which case the volume V2 of zone 2 is larger by a factor 
« 3.3 than the zone 1 volume Vi « 0.05 Vsnr- In a more 
elaborate model, the region between the contact discon- 
tinuity and the reverse shock at 6* « 183" would be inter- 
preted as another zone with a different set of source pa- 
rameters. However, because the prime goal of this paper 
is to establish the need for multi-zone modeling, which is 
appropriate for spatially inhomogeneous sources and can 
significantly relax the constraints on the synchrotron and 
Compton fluxes compared to the single-zone model, here 



we limit our calculations within the framework of two 
zones only. This assumes that zone 2 includes most of the 
shell between the thin X-ray rim (the zone 1) and the re- 
verse shock. This is qualitatively justified as the contact 
discontinuity is apparently porous, with both thermal gas 
and relativistic particles penetrating through the contact 
discontinuity further down, e.g., through Rayleigh- Taylor 
instabilities as evident from detection of thermal X-ray 
filaments [e.g.,[^, 0]. The implied volume filling factors 
of the two zones in the remnant are Ci = ^i/Vsnr ~ 0.05 
and C2 — V2/VSNR ^ 0.5. These values can accomodate 
both interpretations of the X-ray stripes in Tycho [l^l as 
part of either zone 1 (as further inhomogeneities in the 
rim [iOl) or zone 2 (inhomogeneities in the shell). 

The system of coupled equations describing the energy 
distribution functions Ni{E,t) and N2{E,t) of electrons 
(or protons) in a two-zone model was derived in (l^ : 
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Here Pi ee Pi{E,t) and P2 = P2{E,t) represent the 
energy-loss rates {—dE/dt) of particles with energy E 
in zones 1 and 2 respectively. The second term in the 
right side describes convective escape of electrons from 
zone 1 into zone 2, and the 3'"'' term describes the dif- 
fusive escape of electrons from the zone 1 into zone 2 
and vice versa proportional to the difference of the spa- 
tial densities {the gradient) of their energy distributions 
in the transition region between the zones. This leads to 
vanishing of the diffusive exchange of particles between 
the zones at very high energies when the energy densities 
in the two zones become equal. The timescale Tdi{{E) 
depends primarily on the diffusion rate in zone 1 be- 
cause diffusion their is much slower (see fl2| for details). 
Qi = Qi{E,t) and Q2 = Q2{E,t) are the acceleration 
rates of particles in zones 1 and 2, respectively. 

The integral form of the solution for this type of equa- 
tionsegions of can be used to find a self-consistent so- 
lution for the system using the numerical method of it- 
erations [13, neglecting the diffusive influx of particles 
from zone 2 to zone 1 in the first iteration. Calculations 
show that ^ 50 iterations result in a stable (converged) 
solution for Ni{E,t) and N2{E,t). 



III. MODEL PARAMETERS AND RADIATION 
FLUXES FOR TYCHO SNR 

Fluxes of 7-rays produced by relativistic electron 
bremsstrahlung or by hadronic pp-interactions are pro- 
portional to the target gas density Up (in terms of hydro- 
gen atoms) in the source. The gas density ni in zone 1 
(the rim) is proportional to the ambient gas density uq 
upstream of the forward shock, and is about w 4 tiq 
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for a strong shock. Modeling of the X-ray emission 
hnes suggests the best fit ambient gas density po — 
mpUQ ~ 2 X 10~^'*gcm~^, or uq ?s 1.2cm~'^ (see also 
[21]). Smaller values, in the range uq ~ (0.3 — 1) cm^"^ 
for distances d ~ (3 — 2.35) kpc, respectively, have been 
inferred from hydrodynamic modeling of the X-ray emis- 
sion of forward-shocked material [22| . 

The mean gas density n2 in zone 2 (the shell) can 
be derived from the estimate of the total gas mass ac- 
cumulated in zone 2, M2 ~ mpn2C2 Vsnr, where C2 is 
the volume filling factor. Note that in general this 
mass is contributed not only from no, but also from the 
mass of the pre-supernova wind swept up by the for- 
ward shock, and the mass Mq of the SN ejecta. To 
calculate M2 we note that the measured speed of the 
forward shock is significantly smaller than the initial 
vo « IQ^ {Es^.zif/'^ M'^J^kms-^ . For a typical Type 
la SN ejecta with ejecta mass Mq.© = {Mq/Mq) « 1 and 
kinetic energy i^sN.si = ^^SN/lO^^erg « 1 this imphes up 
to several Solar masses residing in the shell of Tycho. 

Indeed, most of the initial explosion energy should 
still be in kinetic form as both the total irradiated en- 
ergy and the energy accumulated in accelerated parti- 
cles are much smaller than i?sN- The total energy of 
the thermal electrons with temperature kTc ~ 1 keV de- 
duced from X-ray observations I, i, is only Et.c = 
{3kTj2)M2/mp ?» 2.5 x 1O4*(M2/M0) erg. To con- 
tribute substantially to the total energy budget in the 
shell, thermal proton plasma should have a temperature 
of order Tp > 10^ K. This is unlikely given fast heat- 
ing of the electrons through Coulomb interactions with 
protons in such two-temperature plasma. Thus, the ki- 
netic energy of the fluid moving with Vf ~ 3wsh/4 behind 
the forward shock in the shell should be still close to the 
initial energy, E^^^ M2vf /2 ~ 10^^i?5i erg. For Vsh 
~ 4600 ((i/2.3) kpc 5] the implied gas density in zone 2 

712 ~ 1.05 (d/3 kpc)-^E5i cm"^ . (3) 

For (2 = 0.5 this leads to 712 in the range from 5 cm^"^ for 
d = 2.5 kpc to ~ 0.9 cm~^ for d = 3.5 kpc. In calculations 
here we use 712 = 3cm~'^ corresponding to 0? « 2.8 kpc. 

For the Compton radiation we take into account not 
only the CMB target photons, but also the FIR photons 
produced by thermal dust in the NE and NW parts of 
Tycho 6]. This component can contribute up to 30% to 
the total Compton flux at TeV energies. 

The particle injection rate Qi is in general time- 
dependent, reaching the maximum at the transition 
time to Sedov phase, and gradually declining afterwards. 
However, taking into account that Tycho is still in its 
early Sedov phase, we assume a stationary injection rate 
Qi{E,t) oc E-^e-^/^'^"' with a = 2.3 and £;cut = 
36 TeV. Given that the reverse shock in Tycho is not 
a bright non-thermal X-ray source, injection of electrons 
in zone 2 at the reverse shock is presumably much smaller 
than at the forward shock. To keep the model simple, we 
neglect electron acceleration in zone 2, i.e. Q2 = 0. 
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FIG. 1: The synchrotron fluxes from radio through X-rays in 
the two-zone model. Dashed and dot-dashed lines show the 
fluxes from zone 1 and zone 2, respectively, and the total flux 
is shown by the solid line. Calculations assume density ri-i ~ 
3cm~"^ at dkpc = 2.8, and n\ ~ n2. Other model parameters 
are Bi = 100 and B2 = 32 ^G, rj = 0.3, a = 2.3 and 
Ecut ~ 36 TeV. Also shown are < GeV bremsstrahlung fluxes 
produced by relativistic electrons in zones 1 and 2. 
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FIG. 2: The 7-ray fluxes produced in the two-zone model. 
The heavy solid line shows the total flux of leptonic origin. 
The total bremsstrahlung and Compton radiation fluxes are 
shown by dashed and solid (thin) lines, respectively. Both 
are predominantly produced in zone 2. For comparison, the 
Compton flux contribution from zone 1 is also shown (dot- 
dashed line). The full open dots show a possible fit to data 
by a hadronic model (see text). 



Figure 1 shows the synchrotron fluxes from radio 
through X-rays calculated assuming Bi = 100 /iC and 
B2 = 32 /iC in zone 1 and zone 2, respectively. In the ra- 
dio to sub-mni wavelengths the total flux is contributed 
mostly by zone 2, while in non-thermal X-rays zone 1 
contributes more than zone 2. This picture is in qualita- 
tive agreement with observations that the narrow rim is 
much brighter in X-rays than in radio. 

Figure 2 shows the 7-ray fluxes. The total flux of 
leptonic origin is dominated by bremsstrahlung up to 
^ 100 GeV, most of which is produced in zone 2 where 
most of electrons accelerated in zone 1 eventually reside. 
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The assumed mean gas density 712 w 3 cm~'^ in zone 
2, at d = 2.8 kpc, implies about 3 AIq of gas accumu- 
lated in the shell of Tycho. The gas density ni in zone 
1 doesn't have any significant impact on the calculated 
fluxes. A density ni in the rim about the same as the 
average n2 in the shell would imply an ambient gas den- 
sity no ~ rii/A = 0.75 cm^'^. Even lower gas densities hq 
cannot, however, be excluded taking into account that 
the mass currently accumulated in the shell may have a 
composition dominated not by the swept up ambient gas 
but by the presupernova wind. 

The total energy of the electrons in zones 1 and 2 are 
Ee.i = 4.4 X 10'''^ erg and £^£.2 = 4.3 x lO'^^erg, respec- 
tively. The thin solid line shows the total Compton flux. 
Because the target photon field (CMBR -I- FIR) pene- 
trates easily through the entire source and thus is same 
in both zones, the Compton contribution from zone I is 
insignificant. 

The detected 7-ray fluxes can also be explained by 
rclativistic protons. Full open dots in Figure 2 show 
the flux from 7r° decay 7-rays produced by protons with 
the power-law index ap = 2.3 and the total energy 
Ep = 3 X 10"*^ erg. We emphasize that the essential dif- 
ference between the hadronic and leptonic fluxes occurs 
only at energies below w 500 MeV. Therefore only by 
detecting the characteristic tt'^ decay cutoff in the energy 
spectrum of 7-rays at £^ ~ (100 — 300) MeV will it be 
possible to claim that cosmic ray hadrons are the origin 
of the 7-rays from Tycho. 

IV. SUMMARY AND DISCUSSION 

Clear identification of the tt" decay bump peaking at 
E — m^o/2 in a photon spectrum was proposed long 
ago as a promising way to establish acceleration sites of 
cosmic-ray protons and nuclei [2^ H^l . The Large Area 
Telescope on Fermi provides a new database for tackling 
this question, and the Type la SNR Tycho presents one of 
the best opportunities to search for signatures of hadronic 
production. 

Any realistic treatment of the acceleration and radia- 
tion must consider at least two zones of different mag- 
netic field strength. Electrons in the thin acceleration 
zone behind the shock where the field is enhanced are 



bright in synchrotron X-rays. Electrons in the region 
with the weaker field (zone 2) occupying the bulk of the 
shell volume will have dim synchrotron emission, and 
can be essentially missed in single-zone analyses. This 
zone can, however, provide most of the Compton and 
bremsstrahlung flux. This effect significantly weakens 
limitations inferred from a single-zone analysis of the 
observed synchrotron X-ray and GeV-TeV 7-ray fluxes. 
With a more realistic two-zone model, we find that a 
good fit to the Fermi and VERITAS data from the Ty- 
cho SNR can be made by a purely leptonic model, with 
bremsstrahlung emission making, primarily, the GeV 7 
rays, and TeV photons coming mainly from Compton- 
scattered soft radiation fields. No hadrons are required 
in this model. The underlying densities and volume filling 
factors of the two zones are not uniquely parameterized, 
our solution represents only one of a range of possibilities. 

Figure 2 shows that hadronic models can fit the data 
just as well as a leptonic model. The two models are 
quite different only at low energies, with the hadronic 
model exhibiting the Tr'^-decay feature, shifted up to a 
few hundred MeV in a lyF^, representation (e.g. psj). 
The absence of this this low energy feature would rule 
out a hadronic origin, while its detection would support 
an origin from cosmic-ray protons and ions. 

Establishing the flux of Tycho at low energies with the 
Fermi-LAT is difficult because of the reduction in the ef- 
fective area and the larger point spread function below 1 
GeV [26*1. Moreover, background modeling of the diffuse 
7-ray emission from the Galaxy disk has to be carefully 
treated. Progress in low-energy analysis and longer expo- 
sure times will help improve the quality of the Fermi-LAT 
spectrum of Tycho. The kinematic tt" flux reduction is, 
however, a generic feature of cosmic-ray hadron origin 
from the GeV cosmic rays that carry most of the energy. 
Identification of this feature in Tycho or other SNRs will 
finally provide experimental confirmation of the SNR ori- 
gin of the Galactic cosmic rays. 
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